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<D . Abstract 

Inclusive-jet and dijet differential cross sections have been measured in neu- 
tral current deep inelastic ep scattering for exchanged boson virtualities Q 2 > 
125 GeV 2 with the ZEUS detector at HERA using an integrated luminosity of 
82 pb _1 . Jets were identified in the Breit frame using the kj< cluster algorithm. 
Jet cross sections are presented as functions of several kinematic and jet vari- 
ables. The results are also presented in different regions of Q 2 . Next-to-leading- 
order QCD calculations describe the measurements well. Regions of phase space 
where the theoretical uncertainties are small have been identified. Measurements 
in these regions have the potential to constrain the gluon density in the proton 
when used as inputs to global fits of the proton parton distribution functions. 
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1 Introduction 



A well established testing ground for perturbative QCD (pQCD) is jet production in 
neutral current (NC) deep inelastic ep scattering (DIS) at high Q 2 , where Q 2 is the 
negative of the square of the virtuality of the exchanged boson V* (V = 7, Z). The 
parton-model process (V*q — > q) in DIS gives rise to final states containing one jet of 
high transverse energy, the so-called 'current jet'. At leading order (LO) in at a , the boson- 
gluon-fusion (V*g — > qq) and QCD-Compton-scattering (V*q — »■ qg) processes give rise to 
two hard jets. 

The hadronic final state in NC DIS may consist of jets with high transverse energy, E^, 
produced in the hard-scattering process, accompanied by the remnant (beam jet) of the 
incoming proton. For the analysis of these types of process, the Breit frame [1] is ad- 
vantageous, since it provides a maximal separation between the products of the beam 
fragmentation and the hard jets. In this frame, the exchanged virtual boson is purely 
space-like, with 3-momentum q = (0,0, —Q), and is collinear with the incoming par- 
ton. In the parton-model process, the virtual boson V* is absorbed by the struck quark, 
which is back-scattered with zero transverse momentum with respect to the V* direction, 
whereas the beam jet follows the direction of the initial struck quark. QCD-Compton 
and boson-gluon-fusion events, on the other hand, are characterised by two partons in 
the final state with possibly non-vanishing transverse momentum with respect to the V* 
direction. Thus, while retaining hard QCD processes at leading order in a s , the contribu- 
tion from the parton-model process can be suppressed by requiring the production of jets 
with high transverse energy in the Breit frame. Therefore, measurements of inclusive- 
jet events at high E^ or dijet events in the Breit frame are directly sensitive to hard 
QCD processes, allowing direct tests of the pQCD predictions and of the proton parton 
distribution functions (PDFs). 

Jet cross sections in NC DIS have been measured previously at HERA. Inclusive-jet [2-4], 
dijet [5] and multijet [6] production have been used to extract values of a s . A previous 
dijet analysis [7] has been used to test the gluon density extracted from global fits. Results 
from inclusive-jet and dijet measurements [3] have been used to constrain the gluon density 
in the proton. Recently, jet cross sections in NC DIS [2] and photoproduction [8] have 
been included in a new NLO QCD analysis to extract the proton PDFs [9], resulting 
in a significant reduction of the uncertainty on the gluon density at medium and high 
momentum fractions. 

This paper presents new measurements of differential dijet cross sections as functions of 
Q 2 , of the mean jet transverse energy of the dijet system in the Breit frame, E?, of the 
dijet invariant mass, Mjj, of the half- difference of the jet pseudorapidities in the Breit 
frame, rj' = {rf^ 1 — f]^ t2 \/2, and of the fraction of the proton momentum taken by the 
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interacting parton, £ = xbj(1 + M 2 JQ 2 ). Here, ^Bj is the Bjorken scaling variable that 
defines, for the parton-model process, the fraction of proton momentum carried by the 
struck parton. Measurements of the dijet cross section as a function of £ are also shown 
for different regions of Q 2 . Also presented are measurements of the inclusive-jet cross 
section as a function of the jet transverse energy in the Breit frame, £^ B! in different 
regions of Q 2 . The analyses are based on data samples with more than twice the statistics 
of the previous studies [2,7]. These measurements probe an extended kinematic regime 
with respect to previous analyses due to the increase in the proton beam energy. The 
improvement in the experimental uncertainties obtained here will facilitate a more precise 
determination of the gluon density in the proton at high £. 



2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [10, 11]. A brief 
outline of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [12], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift-chamber layers, organized in nine superlayers covering the 
polar-angle 1 region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks can be parameterised as (t(j)t)/pt = 0.0058pr © 0.0065 © 0.0014/pr, with px in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the 
energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [13] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and the 
rear (RCAL) calorimeters. Each part is segmented into towers and subdivided in depth 
into one electromagnetic section and either one (in RCAL) or two (in FCAL and BCAL) 
hadronic sections. The smallest subdivision of the calorimeter is called a cell. Under 
test-beam conditions, the CAL single-particle relative energy resolutions were a(E)/E = 
0.18 /^/E for electrons and a(E)/E = 0.35/^/E for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e^p. The 
resulting small-angle energetic photons were measured by the luminosity monitor [14], a 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. The polar angle. 6, and 
the azimuthal angle, (p, are defined with respect to the proton beam direction. The pseudorapidity is 
defined as rj — — In (tan |) . 
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lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

3 Data selection and jet search 

The data used in this analysis were collected during the period 1998-2000, when HERA 
operated with protons of energy E p = 920 GeV and electrons or positrons 2 of energy 
E e = 27.5 GeV, and correspond to an integrated luminosity of 81.7 ± 1.8 pb" 1 . 

A three-level trigger system was used to select events online [11,15]. At the third level, 
NC DIS events were accepted on the basis of the identification of a scattered-electron 
candidate using localised energy depositions in the CAL. An independent trigger selection 
which required at least one (two) jet(s) with transverse energies above 10(6) GeV and 
pseudorapidities below 2.5 was used to measure the fraction of scattered electrons that 
gave a trigger signal. The efficiency of the trigger selection based on the scattered-electron 
candidate was found to be above 99%. 

Events were selected offline using criteria similar to those reported previously [2]. The 
main steps are briefly listed below. The inclusive-jet measurement is based on CAL cells 
for the reconstruction of the kinematic and jet variables, whereas the dijet analysis uses 
a combination of track and CAL information [16]. The selected combinations of tracks 
and CAL clusters are referred to as Energy Flow Objects (EFOs). 

The scattered-electron candidate was identified from the pattern of energy deposits in the 
CAL [17]. The energy (E' e ) and polar angle (8 e ) of the electron candidate were determined 
from the CAL measurements. The kinematic variables Q 2 and xbj were reconstructed 
using the double angle (DA) method [18]. This method uses 9 e and the angle 7^, which 
is equivalent to the angle of the scattered quark in the quark-parton model. The angle 
jh, defined as 

(1 - y)x B jE p - yE e 

COS7h = T\ N 771 , 771 1 

(1 - y)XBiE p + yE e 

was reconstructed using the hadronic final state [18], where y = Q 2 /xBjS and s is the 
centre-of-mass energy. 

The following requirements were imposed on the data sample for both the inclusive-jet 
and dijet analyses: 

• an electron candidate of energy E' e > 10 GeV. This requirement ensured a high and 
well understood electron-finding efficiency and suppressed background from photopro- 
duction events, in which the scattered electron escaped down the rear beampipe; 

2 Here and in the following, the term "electron" denotes generically both the electron and the positron. 
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• the total energy not associated with the electron candidate within a cone of radius 0.7 
units in the pseudorapidity-azimuth (rj — (p) plane around the electron direction should 
be less than 10% of the electron energy. This condition removed photoproduction and 
DIS events in which part of a jet was falsely identified as the scattered electron; 

• the vertex position along the beam axis should be in the range \Z\ < 34 cm. This 
condition helped to select events consistent with ep interactions; 

• -P^miss/ V Et < 2.5 GeV 1 / 2 , where Pr,miss is the missing transverse momentum as 
measured with the CAL and Et is the total transverse energy in the CAL. This cut 
removed cosmic-ray events and beam-related background; 

• Q 2 > 125 GeV 2 ; 

• —0.65 < cos -fh < 0.65. The lower limit avoided a region with limited acceptance due 
to the requirement on the energy of the scattered electron, while the upper limit was 
chosen to ensure good reconstruction of the jets in the Breit frame. 

In addition, the following selection requirements were imposed for the inclusive-jet anal- 
ysis: 

• He < 0.95, where y e = 1 — E' e (l — cos6 e )/(2E e ). This condition removed events in 
which fake electron candidates were found in the FCAL; 

• if the polar angle of the scattered electron was in the range 30° < 9 e < 140°, it was 
required that the fraction of the electron energy within a cone of radius 0.3 units in the 
f] — tp plane around the electron direction should be larger than 0.9; for 9 e < 30°, the 
cut was raised to 0.98. This condition removed events in which jets were misidentified 
as electrons; 

• the event was not consistent with elastic Compton scattering (ep — > ejp), namely no 
second electromagnetic energy cluster above 10 GeV was allowed when the rest of the 
CAL energy, besides the two electromagnetic energy clusters, was below 4 GeV; 

• 38 < (E — Pz) < 65 GeV, where E is the total energy, E = ^ E i: and Pz is the 
Z-component of the vector P = XI? Pi- The sums run over all final-state objects. 
This cut removed events with large initial-state radiation and further reduced the 
background from photoproduction. 

For the dijet analysis, the following conditions had to be fulfilled: 

• Q 2 < 5000 GeV 2 ; 

• 45 < (E - P z ) < 62 GeV. 

The kr cluster algorithm [19] was used in the longitudinally invariant inclusive mode [20] 
to reconstruct jets in the hadronic final state both in data and in Monte Carlo (MC) 
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simulated events (see Section HJ) assuming massless objects. In data, the algorithm was 
applied to the final-state objects after excluding the scattered-electron candidate. The 
jet search was performed in the t] — ip plane of the Breit frame. The jet variables were 
defined according to the Snowmass convention [21]. 

After reconstructing the jet variables in the Breit frame, the massless four-momenta were 
boosted into the laboratory frame, where the transverse energy (-E^lab) an d the pseu- 
dorapidity {rf^ AB ) °f eacn were calculated. Energy corrections [2,8,22,23] were then 
applied to the jets in the laboratory frame and propagated into E^ B in order to com- 
pensate for energy losses in the inactive material in front of the CAL. The following cuts 
were applied: 

• events were removed from the sample if the distance A of any of the jets to the 
electron candidate in the rj — tp plane of the laboratory frame was smaller than 1 unit, 

A = \J (?7lab ~ V e ) 2 + (Vlab ~ V 6 ) 2 < 1> where y? e and if are the azimuthal angle and 
pseudorapidity of the scattered electron, respectively. This requirement removed some 
background from photoproduction and improved the purity of the sample; 

• events were removed from the sample if a jet was in the backward region of the detector 
(?7lab < —2). This requirement removed events in which a radiated photon from the 
electron was misidentified as a jet in the Breit frame; 

• — 2 < rf B < 1.5. This cut restricted the jets to a region with large acceptance. 
The inclusive-jet sample was then selected using the following conditions: 

• -Ey lab > 2.5 GeV. This cut removed a small number of jets for which the uncertainty 
on the energy correction was large; 

• > 8 GeV. 

For the dijet analysis, the following further requirements were imposed: 

• 4 c ; LAB > 3 GeV; 

• Vlab < 2.5. This condition restricted the jets to the region in which the reconstruction 
of the jet variables was optimal; 

• the two highest-Eyg jets in an event, ordered according to E^, were required to 
satisfy E^ > 12 GeV and E^ > 8 GeV. 

The selected sample of inclusive-jet (dijet) events consisted of 19908 (3868) events. 
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4 Monte Carlo simulation and acceptance corrections 



Samples of MC events were generated to determine the response of the detector to jets 
of hadrons and the correction factors necessary to obtain the hadron-level jet cross sec- 
tions. The hadron level is defined in terms of hadrons with lifetime r > 10 ps. The 
generated events were passed through the Geant 3.13-based [24] ZEUS detector- and 
trigger-simulation programs [11]. They were reconstructed and analysed by the same 
program chain as the data. 

Neutral current DIS events including radiative effects were simulated using the Her- 
acles 4.6.1 [25] program with the Djangoh 1.1 [26] interface to the hadronisation 
programs. Heracles includes corrections for initial- and final-state radiation, vertex 
and propagator terms, and two-boson exchange. The QCD cascade is simulated using 
the colour-dipole model (CDM) [27] including the LO QCD diagrams as implemented in 
Ariadne 4.08 [28] and, alternatively, with the MEPS model of Lepto 6.5 [29]. The 
CTEQ5D [30] proton PDFs were used for these simulations. Fragmentation into hadrons 
is performed using the Lund string model [31] as implemented in Jetset 7.4 [32]. 

The jet search was performed on the MC events using the energy measured in the CAL 
cells or EFOs in the same way as for the data. The same jet algorithm was also applied 
to the final-state particles (hadron level) and to the partons available after the parton 
shower (parton level). 

The data were corrected to the hadron level and for QED-radiative effects using bin-by- 
bin correction factors obtained from the MC samples. For this approach to be valid, the 
uncorrected distributions of the data must be well described by the MC simulations. This 
condition was in general satisfied by both the Ariadne and Lepto MC. The Ariadne 
model gave a slightly better description of the inclusive-jet data and was thus used as the 
default model; Lepto was then used to estimate the systematic effect on the correction 
procedure due to the parton-shower model (see Section El). In contrast to this, the dijet 
distributions, shown in Fig. were slightly better described by the Lepto model; here, 
Ariadne was therefore used for systematic checks. In all cases, the correction factors 
differed from unity by typically 10%. These correction factors took into account the 
efficiency of the trigger, the selection criteria and the purity and efficiency of the jet 
reconstruction. 

5 NLO QCD calculations 

The measurements were compared with NLO QCD (0{a 2 s )) calculations obtained using 
the program Disent [33]. The calculations were performed in the MS renormalisation and 
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factorisation schemes using a generalised version [33] of the subtraction method [34]. The 
number of flavours was set to five and the factorisation scale was chosen to be /xp = Q. 
Calculations with different choices of the renormalisation scale, /zr, were performed: the 
default choice was fi R = (E^) 2 for the inclusive-jet analysis; cross checks were performed 
using /j, R = Q 2 and n 2 R = Q 2 + E TB , where E TB is the mean transverse energy of the 

2 

selected jets in an event. For the dijet analysis fj, R = Q 2 + E T was used as default, 
since the resulting predictions describe the data better than either scale alone; here the 
effects of using either Q 2 or E T were investigated. The effects on the NLO calculations of 
using different choices for fi R were found to be smaller than the uncertainties described 
below. The strong coupling constant was calculated at two loops with Al± = 226 MeV, 
corresponding to a s (M z ) = 0.118. The calculations were performed using the CTEQ6 [35] 
parameterisations of the proton PDFs. The kx cluster algorithm was also applied to 
the partons in the events generated by Disent in order to obtain the jet cross-section 
predictions. 

Since the measurements refer to jets of hadrons, whereas the NLO QCD calculations refer 
to jets of partons, the predictions were corrected to the hadron level using the MC models. 
The multiplicative correction factor (Chad) was defined as the ratio of the cross section 
for jets of hadrons over that for jets of partons, estimated by using the MC programs 
described in Section 0] The ratios obtained with Ariadne and Lepto were averaged to 
obtain the value of Chad- This value differs from unity by approximately 5 (5 — 10)% for 
the inclusive-jet (dijet) analysis. 

The NLO QCD predictions for the inclusive-jet analysis were also corrected for the Z°- 
exchange contribution by using MC simulated events with and without Z° exchange. This 
correction is not required for the dijet analysis owing to the restricted Q 2 range. 

Several sources of uncertainty in the theoretical predictions were considered: 

• the uncertainty on the NLO QCD calculations due to terms beyond NLO, estimated 
by varying /i# by a factor of two up and down, was below ±7% at low Q 2 and low E^ B 
and decreased to below ±5% in the highest-Q 2 region for the inclusive-jet analysis. 
For the dijet measurement, the average uncertainty was about ±10%, with the largest 
values approaching ±20% at low £ and low Q 2 values. 

• the uncertainty on the NLO QCD calculations due to that on a s (Mz) was estimated 
by repeating the calculations using two additional sets of proton PDFs, CTEQ6A114 
and CTEQ6A122 [36], determined assuming a s (Mz) = 0.114 and 0.122, respectively. 
The difference between the calculations using these sets and CTEQ6 was scaled by a 
factor of 0.68 to reflect the uncertainty on a s [37]. The resulting uncertainty in the 
cross sections was typically below ±4%; 

• the uncertainty on the modelling of the parton shower was estimated as half the 
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difference between the multiplicative correction factors calculated from the Lepto 
and Ariadne models. The resulting uncertainty on the cross sections was typically 
less than 3%; 

• the uncertainty on the NLO calculations due to the proton PDFs was estimated by 
repeating the calculations using 40 additional sets from the CTEQ6 analysis, which 
takes into account the statistical and correlated systematic experimental uncertainties 
of each data set used in the determination of the proton PDFs. The resulting uncer- 
tainty in the cross sections was about ±4% at low Q 2 and decreased to around ±2% 
at high Q 2 ; 

• the uncertainty of the calculations in the value of fj,p was estimated by repeating the 
calculations with /xp = Q/2 and 2Q. The variation of the calculations was negligible. 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. Figures 121 and El show an overview of various theoretical uncer- 
tainties for both the inclusive-jet and the dijet analysis. 

The gluon-induced contribution to the cross section as a function of ^ et B for the inclusive- 
jet analysis and of log 10 £ for the dijet analysis in various regions of Q 2 is shown in Figs. |U 
and 03 The gluon fraction varies between 70% at low Q 2 for both low £ and low E^ B 
and less than 10% for the lowest E^ B at high values of Q 2 . For the dijet analysis, the 
gluon contribution is always at least 30%. In addition, Fig. shows, separately for the 
inclusive-jet and dijet analyses, the gluon-induced contribution to the cross section as a 
function of Q 2 together with the uncertainties as taken from the CTEQ6 analysis. 

6 Experimental uncertainties 

The following sources of systematic uncertainty were considered for the measured cross 
sections: 

• the uncertainty in the absolute energy scale of the jets was estimated to be ±1% for 
■^t lab > 10 GeV and ±3% for lower E^ LAB values [8,22,38]. The resulting uncertainty 
on the cross sections was about ±5% and increased to approximately ±10% in certain 
regions of the dijet phase space; 

• the uncertainty in the absolute energy scale of the electron candidate was estimated 
to be ±1% [39]. The resulting uncertainty was below ±1%; 

• the differences in the results obtained by using either Ariadne or Lepto to correct 
the data for detector effects were typically below ±8%; 

• the analysis was repeated using an alternative technique [40] to select the scattered- 
electron candidate. The resulting uncertainty was typically below ±3%; 
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• the -E^lab cut was raised to 4(4.5) GeV for the inclusive-jet (dijet) analysis. The 
resulting uncertainty was typically smaller than ±1%; 

• the cut in rf^ AB used to suppress the contamination due to photons falsely identified 
as jets in the Breit frame was changed to —3 and to —1.5. The resulting uncertainty 
was typically below ±1%; 

• the uncertainty due to the selection cuts was estimated by varying the values of the cuts 
within the resolution of each variable. The effect on the cross sections was typically 
below ±2%; 

• the uncertainty on the cross sections due to that in the simulation of the trigger was 
negligible. 

The systematic uncertainties not associated with the absolute energy scale of the jets were 
added in quadrature to the statistical uncertainties and are shown in figures l7ltolT2l as error 
bars. The uncertainty due to the absolute energy scale of the jets is shown separately as a 
shaded band in each of these figures, due to the large bin-to-bin correlation. In addition, 
there was an overall normalisation uncertainty of 2.2% from the luminosity determination, 
which is not included in the figures. 

7 Results 

7.1 Dijet differential cross sections 

Dijet cross sections were measured in the kinematic region 125 < Q 2 < 5000 GeV 2 and 
| cos 7^ | < 0.65. These cross sections correspond to the two highest-£ , ^ t B jets in each event 
with E$tb > 12 GeV > 41 > 8 GeV and - 2 < # < L5 and 

were corrected for detector 

and QED radiative effects as described in Section 0] 

The measured cross sections as functions of Q 2 , x-Qj, E t , Mjj, rj' and log 10 £ are shown 
in Figs. and |H1 and are listed in Tables Q] to El In these and subsequent figures, each 
data point is plotted at the abscissa at which the differential cross section was equal to 
its bin-averaged value, according to the NLO QCD calculation. The data distribution 
as a function of Q 2 exhibits a fall-off of more than two orders of magnitude within the 
range studied. The Q 2 -restricted range implicitly limits the £ range of the measurements 
to 0.009 < £ < 0.37, a region of phase space where these measurements can significantly 
constrain the gluon density in the proton. The dijet cross sections as functions of Et and 
My are particularly suited to test the matrix elements in the perturbative calculations. 
The measured distributions show a steep fall-off of more than two orders of magnitude 
within the measured range. 
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The NLO QCD calculations with n 2 R = Q 2 + E T are compared to the data in Figs. [7| 
and |H1 They give a good description of the data. The calculations with \j? R = Q 2 or 
fi 2 R = E T are included in Figs. [7| and |S] as well. At low Q 2 , low Mj and low Et, the 
theoretical uncertainty is dominated by the effect of the variation of the renormalisation 
scale. In the case of the dijet cross sections as functions of Et and My v the theoretical 
uncertainty is smaller than the experimental at high values of Et and My. The excellent 
agreement between data and theory for these distributions demonstrates the validity of 
the description of the dynamics of dijet production by pQCD at 0(aa 2 ). 

7.2 Dijet and inclusive-jet differential cross sections in different 
regions of Q 2 

Figure El shows the measured dijet cross section as a function of log 10 £ in different regions 
of Q 2 . The cross sections are also given in Table [7| The requirement that two jets with 
high transverse energy be observed in the final state suppresses the cross section in the 
low-£ region, so that the measured cross section rises as log 10 £ increases. At high values 
of log 10 £, the decrease of the cross section reflects the decrease of the gluon and quark 
densities at high £. 

The NLO QCD predictions with /i R = Q 2 + E T are compared to the measurements in 
Fig-EHusing the relative difference of the measured differential cross sections to the NLO 
calculations; the uncertainty of the calculation is also shown in the figure. The predictions 
of NLO QCD give a good description of the data. Also indicated are the predictions using 
/4 = Q 2 or /i| = E T . 

Inclusive-jet cross sections were measured in the kinematic region Q 2 > 125 GeV 2 and 
ICOS7/J < 0.65. These cross sections include every jet of hadrons in the event with 
£$ B > 8 GeV and -2 < rf B < 1.5 and were corrected for detector and QED radiative 
effects. The cross sections for different regions of Q 2 as a function of E^ B are presented in 
Fig. ^2 and Table |H1 The measured cross sections exhibit a steep fall-off within the E^ B 
range considered. As Q 2 increases, the £^ ot B dependence of the cross section becomes less 
steep. 

The NLO QCD predictions with ji R = (E 3 ^) 2 are compared to the measurements in 
Fig. HH They give a good description of the data. To study the scale dependence, 

2 

calculations using fi 2 R = Q 2 or fJ, 2 R = Q 2 + E T B are also compared to the data in Fig. ITTI 
they provide an approximately equally good description of the data. 

Figure shows the relative difference of the measured differential cross sections to the 
NLO QCD calculations with /i R = (E^) 2 . The uncertainty of the NLO QCD calculations 
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and the results of the theory calculations using either = Q 2 or n 2 R = Q 2 + E TB are 
also shown. The data are well described by the predictions. 

As indicated in Section the theoretical uncertainties in both the inclusive-jet and dijet 
regimes are dominated by the contribution from higher orders as estimated by varying 
the renormalisation scale. As was shown in Figs. 121 and H3 the scale uncertainty decreases 
as Q 2 increases; it is also lower for high values of £ than for low values of this variable. 
The contribution from the PDF uncertainty is approximately constant in all variables and 
non-negligible. Especially at high values of £ (E 3 ^) for the dijet (inclusive-jet) analysis, 
the PDF uncertainty is even dominant. In these regions, in which the gluon-induced 
contribution is still substantial (see Figs. 0] to the data will be able to constrain 
further the gluon density in the proton. 

8 Summary and conclusions 

Inclusive-jet and dijet cross sections have been measured at high Q 2 values at HERA 
for yt~s = 318 GeV. The jets were reconstructed using the kx cluster algorithm in the 
longitudinally invariant inclusive mode in the Breit frame. The dijet cross sections were 
measured as functions of Q 2 , xbj, E t , Mjj, rj' and £. In addition, inclusive-jet and dijet 
measurements were performed in different regions of Q 2 as functions of E^ B and £, respec- 
tively. The data are well described by the NLO QCD predictions. The cross sections in 
different regions of Q 2 are shown to be sensitive to the gluon density of the proton. The 
precise measurements obtained here are therefore of particular relevance for improving 
the determination of the gluon density in future QCD fits. 
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Q 2 bin 
(GeV 2 ) 


da/dQ 2 

(pb/GeV 2 ) 8 Btat 4 y st 5es 


Cqed 


Chad 


125-250 

250-500 

500-1000 

1000-2000 

2000-5000 


0.1183 ±0.0033 ±g : ggS 
0.0589 ±0.0018 IgX IS 
0.02061 ±0.00074 t° Zlt ±Sl8SSS 
0.00602 ±0.00028 t° Zll IE 
0.001189 ±0.000075 t° Z°otl +°oZ°oll 


0.96 
0.94 
0.93 
0.91 
0.96 


0.85 

0.91 
0.92 
0.95 
0.96 



Table 1: Dijet cross-section da/dQ 2 for jets of hadrons in the Breit frame selected 
with the longitudinally invariant kr cluster algorithm. The statistical, uncorrelated 
systematic and jet-energy-scale (ES) uncertainties are shown separately. The mul- 
tiplicative corrections applied to the data to correct for QED radiative effects, Cqed, 
and the corrections for hadronisation effects to be applied to the parton-level NLO 
QCD calculations, Chad? o,re shown in the last two columns. 



xbj bin 


da/dx Bj 

(pb) 5 s tat #syst #ES 


Cqed 


Chad 


0.0001-0.01 
0.01-0.02 
0.02-0.035 
0.035-0.07 
0.07-0.1 


1599 ±45 tit tlf 

1849 ±48 tH -95 

647 ±24 tit -It 
141 7 +7 e +4.0 +2.9 

91 1 10 9 +3.3 +0.5 
Zl.l HZO.Z _ 33 _ Q1 


0.95 
0.94 
0.92 
0.95 
0.89 


0.86 
0.93 
0.92 
0.91 
0.90 



Table 2: Dijet cross-section dojdx^ for jets of hadrons in the Breit frame se- 
lected with the longitudinally invariant kr cluster algorithm. Other details as in 
the caption to Tabled 
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Et bin 


dcr/dET 












(GeV) 


(pb/GeV) 


^stat 


^syst 


<$ES 


Cqed 


Chad 


10-16 


5.71 


±0.11 


+0.12 
-0.12 


+0.30 
-0.26 


0.94 


0.90 


16-22 


1.935 


±0.064 


+0.086 
-0.084 


+0.081 
-0.091 


0.94 


0.90 


22-30 


0.361 


±0.023 


+0.023 
-0.023 


+0.017 
-0.016 


0.94 


0.91 


30-60 


0.0371 


±0.0043 


+0.0076 
-0.0076 


+0.0021 
-0.0018 


0.97 


0.91 



Table 3: Dijet cross-section da/dE? for jets of hadrons in the Breit frame se- 
lected with the longitudinally invariant k? cluster algorithm. Other details as in 
the caption to Tabled 



Mjj bin 
(GeV) 


da/dMjj 

(pb/GeV) <5 stat <5 sy st S ES 


Cqed 


Chad 


20-32 
32-45 
45-65 
65-120 


2 382 ±0 051 + 055 +0 - 123 

Z.OOZ 3ZU.U01 -0.054 -0.109 

1 1 OA +n OQ4 +0-047 +0.049 
1.10+ ICU.UO^ -0.047 -0.055 

n 999 _i_nni9 +0.006 +0.009 

u -^ z zcu.Uiz, _o.006 -0.009 

0.0180 ±0.0021 ±K? + -°oZI 


0.95 
0.94 
0.92 
0.97 


0.91 

0.89 
0.89 
0.94 



Table 4: Dijet cross-section da/dM^ for jets of hadrons in the Breit frame se- 
lected with the longitudinally invariant kr cluster algorithm. Other details as in 
the caption to TableUi 



rj' bin 


da / dr)' 

(pb) 5 s tat 4yst S E S 


Cqed 


Chad 


0-0.1 
0.1-0.25 
0.25-0.45 
0.45-0.65 
0.65-1.60 


87.9 ±3.4 t|° i|? 
87.0 ±2.8 ±\i tf 6 7 
68.5 ±2.1 t\i t|o 
42.7 ±1.6 lt° tH 

6.07 ±0.28 t° % •!;•;! 


0.95 
0.94 
0.95 
0.94 
0.91 


0.91 

0.94 
0.91 
0.87 
0.86 



Table 5: Dijet cross-section da/dr/' for jets of hadrons in the Breit frame se- 
lected with the longitudinally invariant kx cluster algorithm. Other details as in 
the caption to Tabled 
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logio £ bin 


(Pb) 5 sta t ^syst <$ES 


Cqed 


Chad 


-2 - -1.5 
-1.5 - -1.35 
-1.35 - -1.1 
-1.1 - -0.85 
-0.85 - -0.5 


22.22 ±0.78 to 7 7 l t{f 7 
74.4 ±2.6 tH ±11 
73.9 ±1.9 ±« j** 
31.0 ±1.3 t\i tH 

Q QQ +n Qfi +0.48 +0.32 

o.yo zcu.ou _ 04g _ 021 


0.95 
0.94 
0.94 
0.94 
0.93 


0.88 

0.92 
0.92 
0.87 
0.79 



Table 6: Dijet cross-section da/dlog l0 ^ for jets of hadrons in the Breit frame 
selected with the longitudinally invariant kr cluster algorithm. Other details as in 
the caption to Tabled 
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logic 


£ bin 


da/d\og lQ £ 
















(pb) 


Ostat 


^syst 




Cqed 


Chad 


125 < Q 2 < 250 GeV 2 


-2.00 


- -1.50 


9.17 


±0.48 


+0.44 
-0.44 


+0.70 
-0.62 


0.97 


0.83 


-1.50 


- -1.35 


25.0 


±1.4 


+1.8 
-1.8 


+1.4 
-0.7 


0.95 


0.89 


-1.35 


- -1.10 


19.74 


±0.93 


+0.43 
-0.41 


+1.17 
-1.65 


0.95 


0.86 


-1.10 


- -0.50 


2.52 


±0.20 


+0.10 
-0.10 


+0.26 
-0.19 


0.99 


0.83 


250 < Q 2 < 500 GeV 2 


9 nn 
— z.uu 


- — l.OU 


8.65 


±0.49 


+0.59 
-0.57 


+0.42 
-0.38 


n Q/i 


u.yo 


-1.50 


- -1.30 


23.4 


±1.2 


+1.5 
-1.4 


+1.1 
-1.3 


0.94 


0.93 


-1.30 


- -1.00 


16.18 


±0.83 


+0.55 
-0.44 


+1.02 
-0.78 


0.93 


0.89 


-1.00 


- -0.50 


1.77 


±0.21 


+0.09 
-0.09 


+0.10 
-0.11 


0.94 


0.88 






500 < Q 


2 < 1000 GeV 2 








-1.90 


- -1.50 


4.32 


±0.40 


+0.15 
-0.18 


+0.23 
-0.17 


0.94 


0.92 


-1.50 


- -1.20 


17.50 


±0.88 


+1.00 
-1.00 


+0.72 
-0.75 


0.94 


0.94 


-1.20 


- -0.90 


9.45 


±0.63 


+0.90 
-0.90 


+0.50 
-0.35 


0.91 


0.92 


-0.90 


- -0.60 


1.48 


±0.26 


+0.09 
-0.09 


+0.06 
-0.06 


0.98 


0.86 






1000 < Q 2 < 2000 GeV 


2 






-1.70 


- -1.40 


2.72 


±0.37 


+0.33 
-0.33 


+0.06 
-0.12 


0.93 


0.94 


-1.40 


- -1.20 


10.98 


±0.86 


+1.23 
-1.23 


+0.22 
-0.01 


0.94 


0.98 


-1.20 


- -1.00 


10.37 


±0.84 


+0.48 
-0.48 


+0.26 
-0.28 


0.92 


0.96 


-1.00 


- -0.60 


2.17 


±0.25 


+0.15 
-0.15 


+0.06 
-0.08 


0.85 


0.89 






2000 < Q 2 < 5000 GeV 


2 






-1.60 


- -1.20 


2.14 


±0.26 


+0.14 
-0.14 


+0.01 
-0.02 


0.92 


0.99 


-1.20 


- -1.00 


8.28 


±0.80 


+0.54 
-0.54 


+0.14 
-0.18 


1.00 


0.96 


-1.00 


- -0.60 


2.48 


±0.30 


+0.16 
-0.16 


+0.08 
-0.04 


0.93 


0.94 



Table 7: Dijet cross-sections da/d\og l0 l; for jets of hadrons in the Breit frame 
selected with the longitudinally invariant cluster algorithm in different regions 
of Q 2 . Other details as in the caption to Tabled 
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Ej,-q bin 


da/dEj.-Q 












(GeV) 


(pb/GeV) 


"Stat 


^syst 


<5es 


Cqed 


Chad 




125 < Q 2 


< 250 GeV 2 






8-10 


32.97 


0.49 


+1.21 
-1.21 


+1.81 
-1.69 


0.96 


0.90 


10-14 


13.00 


0.22 


+0.19 
-0.19 


+0.79 
-0.75 


0.98 


0.94 


14-18 
18-25 


3.71 
0.835 


0.11 
0.037 


+0.15 
-0.15 
+0.013 
-0.012 


+0.28 
-0.24 
+0.051 
-0.056 


0.97 
0.94 


0.94 
0.93 


25-100 


0.0160 


0.0014 


+0.0027 
-0.0027 


+0.0010 
-0.0011 


0.97 


0.86 




250 < Q 2 


< 500 GeV 2 






8-10 


18.40 


0.38 


+0.74 
-0.74 


+0.68 
-0.60 


0.94 


0.92 


10-14 


8.74 


0.19 


+0.30 
-0.30 


+0.33 
-0.35 


0.96 


0.95 


14-18 
18-25 


3.30 
0.889 


0.11 
0.042 


+0.15 
-0.15 
+0.041 
-0.041 


+0.18 
-0.14 
+0.052 
-0.057 


0.96 
0.92 


0.97 
0.97 


25-100 


0.0242 


0.0020 


+0.0005 
-0.0005 


+0.0012 
-0.0011 


0.95 


0.91 


500 < Q 2 < 1000 GeV 2 


8-10 


8.79 


0.26 


+0.34 
-0.34 


+0.26 
-0.15 


0.96 


0.91 


10-14 
14-18 


4.69 
2.239 


0.14 
0.093 


+0.19 
-0.19 
+0.137 
-0.137 


+0.11 
-0.13 
+0.091 
-0.074 


0.94 
0.93 


0.95 
0.98 


18-25 


0.701 


0.039 


+0.051 
-0.051 


+0.026 
-0.026 


0.96 


0.99 


25-100 


0.0335 


0.0027 


+0.0018 
-0.0018 


+0.0019 
-0.0019 


0.96 


0.97 




1000 < Q 2 


< 2000 GeV 2 






8-10 
10-14 


3.30 
1.985 


0.16 
0.091 


+0.14 
-0.14 
+0.077 
-0.077 


+0.09 
-0.08 
+0.017 
-0.037 


0.93 
0.91 


0.93 
0.95 


14-18 
18-25 


1.115 
0.492 


0.069 
0.034 


+0.056 
-0.056 
+0.039 
-0.039 


+0.034 
-0.008 
+0.009 
-0.018 


0.98 
0.93 


0.99 
0.99 


25-100 


0.0263 


0.0026 


+0.0043 
-0.0043 


+0.0020 
-0.0013 


1.00 


1.00 




2000 < Q 2 


< 5000 GeV 2 






8-10 


1.292 


0.095 


+0.120 
-0.120 


+0.033 
-0.022 


0.92 


0.90 


10-14 


0.858 


0.060 


+0.024 
-0.024 


+0.007 
-0.006 


0.90 


0.93 


14-18 


0.612 


0.052 


+0.070 
-0.070 


+0.023 
-0.017 


1.02 


1.00 


18-25 
25-100 


0.242 
0.0185 


0.024 
0.0021 


+0.028 
-0.028 
+0.0023 
-0.0023 


+0.009 
-0.006 
+0.0004 
-0.0005 


0.96 
0.91 


1.00 
0.99 


5000 < Q 2 < 100000 GeV 2 


8-10 


0.225 


0.037 


+0.091 
-0.091 


+0.011 
-0.006 


0.99 


0.93 


10-14 
14-18 
18-25 


0.267 
0.122 
0.070 


0.037 
0.024 
0.013 


+0.022 
-0.023 
+0.017 
-0.017 
+0.019 
-0.019 


+0.003 
-0.019 
+0.003 
-0.005 
+0.001 
-0.000 


0.96 
0.97 
0.98 


0.93 
0.98 
0.99 


25-100 


0.0114 


0.0022 


+0.0042 
-0.0042 


+0.0003 
-0.0001 


0.99 


1.00 



Table 8: Inclusive jet cross-sections da/dE^ B for jets of hadrons in the Breit 
frame selected with the longitudinally invariant cluster algorithm in different 
regions of Q 2 . Other details as in the caption to Tabled 
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Figure 1: Uncorrected data distributions for dijet production with E^ B > 12 GeV, 
-E^b > 8 GeV and —2 < rf^ < 1.5 in the kinematic range given by \cosjh\ < 0.65 
and 125 < Q 2 < 5000 GeV 2 (dots). For comparison, the predictions of the ARI- 
ADNE (dashed histograms) and LEPTO (solid histograms) MC models are also 
included. 
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Figure 2: Overview of theoretical uncertainties for the inclusive- jet analysis as 
functions of E^ B in bins of Q 2 for Icc^^l < 0.65 and —2 < r/g* < 1.5. Shown are 
the relative uncertainties induced by the variation of the renormalisation scale hr 
(hatched area) and by the uncertainties on the proton PDFs (shaded area). 
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Figure 3: Overview of theoretical uncertainties for the dijet analysis as functions 
o/log 10 £ in bins of Q 2 for |cos7 h | < 0.65, E J T e £ > 12 GeV, E 3 T e g > 8 GeV and 

—2 < T]^* < 1.5. Other details as in the caption to Fig. [B 
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Figure 4: NLO QCD predictions of the gluon-induced fraction of the inclusive- 
jet cross sections da/dE^ as functions of E^ B in different regions of Q 2 for 

IC0S7/J < 0.65 and —2 < rf# < 1.5 (shaded histograms). The CTEQ6 parameteri- 
sations of the proton PDFs were used. 
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Figure 5: NLO QCD predictions of the gluon-induced fraction of the dijet 
cross-sections da/d\og 10 £ as functions o/log 10 £ in different regions of Q 2 for 
\cos lh \ < 0.65, E^l > 12 GeV, > 8 GeV and -2 < r/ B ct < 1.5. Other 

details as in the caption to Fig. ^| 
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Figure 6: NLO QCD predictions of the gluon-induced fraction of the di- 
jet (inclusive-jet) cross-section da/dlog^ (da/dE 3 ^) as a function of Q 2 for 

|cos7 h | < 0.65 and -2 < r/ B ct < 1.5 and for log 10 £ = -1.25 (10 < E^ B < U GeV); 
the uncertainty on the gluon-induced fraction due to those on the proton PDFs is 
represented by the shaded (hatched) area. The CTEQ6 parameterisations of the 
proton PDFs were used. 
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Figure 7: The measured differential cross-sections (a) da/dQ 2 and (b) dajdx^ 
for dijet production with -E^b > ^ GeV, -E^b > & GeV and —2 < rj 3 ^ < 1.5 
(dots), in the kinematic range given by |cos7h| < 0.65 and 125 < Q 2 < 5000 GeV 2 . 
The inner error bars represent the statistical uncertainty. The outer error bars 
show the statistical and systematic uncertainties, not associated with the uncer- 
tainty in the absolute energy scale of the jets, added in quadrature. The shaded 
bands display the uncertainties due to the absolute energy scale of the jets. The 
NLO QCD calculations with [x 2 R = Q 2 + E T (solid lines), jj? R = Q 2 (dashed lines) 

and n 2 R = E T (dash-dotted lines), corrected for hadronisation effects and using the 
CTEQ6 parameterisations of the proton PDFs, are also shown. The lower parts of 
the figures show the relative differences with respect to the NLO QCD calculations 

2 

with (i 2 R = Q 2 + E T : measurements (dots), NLO QCD calculations with (i 2 R = Q 2 

2 

(dashed lines) and with (i 2 R = E T (dotted lines); the hatched bands display the total 
theoretical uncertainty. 
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Figure 8: The measured differential cross-sections (a) da/dE T , (b) da/dM^, (c) 
da/drj' and (d) da/dlog 10 ^ for dijet production with E^ B > 12 GeV, -E^b > & GeV 
and —2 < r/g* < 1.5 (dots), in the kinematic range given by |cos7/j| < 0.65 and 
125 < Q 2 < 5000 GeV 2 . Other details as in the caption to Fig. [7| 
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Figure 9: The measured differential cross-section da/d\og 10 ^ for dijet production 
with -E^b > 12 GeV, -E^b > 8 GeV and —2 < rf^ < 1.5 in different regions of Q 2 
(dots). Other details as in the caption to Fig.^ 
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Figure 10: Relative differences between the measured differential cross- sections 
da/d log 10 £ presented in Fig.\Qand the NLO QCD calculations with n R = Q 2 + E T 
(dots). The relative differences between the NLO QCD calculations with = Q 2 

(n R = E T ) and those with /j iR = Q 2 + E T are also shown as dashed (dotted) lines. 
Other details as in the caption to Fig. ^ 
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Figure 11: The measured differential cross-section da/dE ] f l B for inclusive-jet 
production with > 8 GeV and —2 < rf B < 1.5 in different regions of Q 2 

(dots). The NLO QCD calculations with ff 2 ^ = (B^g) 2 (solid lines), [i 2 R = Q 2 

(dashed lines) and /j, r = Q 2 + E TB (dash-dotted lines), corrected for hadronisation 
and Z° effects and using the CTEQ6 parameterisations of the proton PDFs, are 
also shown. Other details as in the caption to Fig. ^ 
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Figure 12: Relative differences between the measured differential cross-sections 
da/dE^s presented in Fig. \H\ and the NLO QCD calculations with n R = (E^ B ) 2 
(dots). The relative differences between the NLO QCD calculations with n 2 R = Q 2 



(v^r = Q 2 + E TB ) and those with fi 2 R = (Ep B ) 2 are also shown as dashed (dash- 
dotted) lines. Other details as in the caption to Fig. 
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